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+L Abstract 

K/1 I employ the van der Waals theory of Baus and coworkers to analyze the fast, adiabatic decay 

2 of a supercooled liquid in a closed vessel with which the solidification process usually starts. By 

I imposing a further constraint on either the system volume or pressure, I use the maximum-entropy 

o 
o 



method to quantify the fraction of liquid that is transformed into solid as a function of undercooling 
and of the amount of a foreign gas that could possibly be also present in the test tube. Upon 
looking at the implications of thermal and mechanical insulation for the energy cost of forming a 



^1^ solid droplet within the liquid, I identify one situation where the onset of solidification inevitably 



occurs near the wall in contact with the bath. 
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I. INTRODUCTION 

Liquid freezing is a widely studied phenomenon, especially under equilibrium conditions 
where the temperature T and the pressure P of the system are kept fixed through the 
contact with a bath. Considerable attention has also been (and still currently is) devoted to 
the kinetic aspects of freezing, which in a moderately supercooled liquid is initiated by the 
spontaneous nucleation and growth of a sizable crystal droplet. The usual setting where the 
initial stages of freezing are studied is again isothermal-isobaric, which assumes a prompt 
release of the latent heat of solidification to the system environment. In the present study, I 
keep distinct the two steps by which the freezing process actually develops for a supercooled 
liquid [H 12] , namely a rapid return to the solid-liquid equilibrium temperature (this stage 
occurs so quickly that there would be no time for a significant transfer of heat to the bath), 
followed by the slower, diathermal solidification process governed by heat conduction to the 
bath, which, in effective terms, begins in a later moment. While the late freezing stage has 
been the focus of many studies (inspiring a whole branch of mathematical physics which 
goes under the name of "Stefan problem", see e.g. [3|), the initial adiabatic-freezing process 
has received less attention in the literature, being confined to applied-research areas like 
Atmospheric Science (where it is studied in connection with the physics of ice accretion 
and chemicals uptake in hydrometeors |1HZ]) and Metallurgy (bearing here the name of 
"recalescence" , see e.g. [SIE])- 

Clearly, the initial adiabatic transformation to a two-phase state can only be observed 
if the process of solid growth takes a time teq much smaller than the time tdia that would 
be needed to transfer (essentially by conduction) an energy equal to the latent heat of 
complete freezing to the surroundings - differently, solidification would occur without this 
intermediate two-phase state. This can be represented as 

teq^y^/VC^tdia, (1.1) 

where V is the system volume, G is the long-term growth rate of the solid inside the liq- 
uid (measured e.g. in cm/s and likely proportional to the liquid supersaturation for small 
undercoolings), and tdia is the diathermal-freezing time (see Eq. (2) of Ref. pj for an esti- 



mate). Whether Eq. (1.1) is satisfied or not in a concrete case, it would not be easy to say 
(see, however, below concerning the consequences of thermal insulation for the energetics of 
droplet formation and the kinetics of crystallization). 



At the end of the adiabatic stage, a microsegregated solid-liquid mixture with definite pro- 
portions of the two phases is formed, which later undergoes complete crystallization through 
the exchange of energy and possibly volume with the colder bath [TUHT3] . A fine-grained, 
uniform distribution of the solid within the liquid is a frequent outcome of recalescence (see 
e.g. the discussion in [14J), as confirmed experimentally for water by Hindmarsh et al. [2] 
and well known to metallurgists. When a molten alloy (or a multicomponent liquid what- 
soever) is cooled, the growing solid phase usually forms a porous matrix through which the 
residual liquid can fiow. The reactive medium made of the solid matrix and the residual 
liquid is called a mushy zone. This type of semisolid, solid-liquid mixture occurs for many 
freezing conditions forming dendrites, see e.g. [I5]. The full conversion of the mushy zone 
into a compact solid occurs through a coarsening mechanism (also called Ostwald ripening), 
involving mass and heat diffusion, and requires some time during which the temperature of 
the system stays constant at the equilibrium freezing temperature. 

From a thermodynamic point of view, the driving force of the return of the metastable 
liquid to equilibrium is entropy maximization, insofar as the system and (when present) the 
volume reservoir are treated as parts of a larger thermally isolated system. For this system, 
the final equilibrium state is the one with the maximum entropy under the given internal 
constraints |16]. In the mushy zone, however, entropy is not yet at a maximum because of 
the large amount of interfacial energy trapped in the interstices of the solid matrix. Hence, 
the state of the system immediately after recalescence is not the one prescribed by ther- 
modynamics but rather a compromise dictated by kinetic considerations. While not being 
completely realistic, a study of adiabatic freezing solely based on the maximum-entropy 
principle will anyway help to analyze general trends of thermodynamic variables with su- 
percooling. This would be especially true for rare-gas fiuids, having a weakly anisotropic 
solid-liquid interface tension, and in the low-supercooling regime where, due to a low nucle- 
ation rate, the size attained by solid grains before impinging other grains will be larger. 

In maximizing entropy, besides the adiabaticity of the system boundaries, additional 
constraints to be accounted for may regard the total system mass (which is here assumed 
to be conserved) and its volume or pressure, depending on the physical context. To make 
progress, I shall represent the system characteristics by a specific model of simple-fiuid 
thermodynamics, viz. the phenomenological theory of Bans and coworkers [T71 HH] , which, 
though being of little quantitative value, at least provides a reasonable close-form entropy 



function for each phase. With this tool at hand, I shall demonstrate that the equilibrium 
state eventually attained after completion of the adiabatic process is indeed inhomogeneous, 
with both solid and liquid present in calculable proportions. 

This paper is organized as follows. After introducing the model and the method in Section 
II, I analyze three different physical situations of adiabatic freezing in Section III, for each 
deriving a number of numerical results. Then, in Section IV these results are exploited to 
see whether the energy cost associated with the formation of the solid-liquid interface could 
be a stumbling block to crystallization under isolated conditions. Some concluding remarks 
are presented in Section V. 

II. MODEL AND METHOD 

The simplest set-up for the study of adiabatic freezing is the following. Consider a A^- 
particle liquid which completely fills a closed rigid container of volume V. The liquid, initially 
in stable equilibrium at the melting/freezing temperature T^, is then gently cooled until the 
temperature Tin < T^ is reached. At this point, the pressure and energy of the liquid (Pm 
and E[^) can be accessed, at least in principle, from the metastable branches of its mechanical 
and thermal equations of state. Now imagine to remove the contact with the thermostat 
and to induce then (e.g. by a mechanical shock) the irreversible decay of the system to 
equilibrium (this transformation is both adiabatic and isochoric, hence energy-conserving). 
We want to determine in which equilibrium state the system will eventually settle down. 
This question may be answered by appealing to the maximum-entropy principle. Envisaging 
the possibility that part of the system could remain liquid, I denote by Ei,Vi, and Ni the 
energy, volume, and particle number of the liquid fraction. If the fundamental relations 
(entropy functions) of the liquid and solid phases are known, the total system entropy reads 

Stot = Si{Ei, Vi, Ni) + 5,(Ein -Ei,V-Vi,N- Ni) , (2.1) 

assuming weak coupling between the phases and taking the internal constraints into account. 
The necessary conditions for the maximum of Stot are three, one for each liquid state vari- 
able, and are equivalent to requiring that the temperature, the pressure, and the chemical 
potential of the solid and liquid components be equal. We then see that (unless Ni = 



or A^ at the point of maximum of (2.1)) the final equilibrium state lies on the solid-liquid 



coexistence locus, as empirically observed. In order for Eq. (2.1 ) to be really useful, however, 
the functions Si and 5*^ are to be made explicit. 

The mean-field theory by Bans and coworkers [T71 [12] provides a convenient framework 
for discussing adiabatic freezing at a semiquantitative level. It is meant for a system of 
particles interacting through a spherically-symmetric potential u{r) given by the hard-sphere 
potential plus a short-range attractive tail, 

{+00 , r < a 

. . ^ ' ^2.2) 

—t(j)[r /a) , r > a 

where a is the particle-core diameter, e is the depth of the attractive well, and 0(x) > 
gives the shape of the well in terms of the scaled interparticle distance x = rjo. In the same 
spirit of the van der Waals theory, the repulsive and attractive potential terms separately 
concur to build up the system Helmholtz free energy, which is taken to be 

F(T, V, N) = Fn{T, V, N) + Fa{T, V, N) , (2.3) 

where -Fr(T, V, N) = Fid(T, aV, N) is the free energy of N fictitious non-interacting particles 

in a fraction a of the total system volume, and 

AT r 
Fa{T,V,N) = - j dVp(r)w^(r), (2.4) 

where p(r) is the local number density experienced by a reference particle in the origin and 
ua is the attractive potential. The description of the model is complete after specifying a 
and p(r) for each phase. For a fluid phase, one assumes 

p , _, . N 

a = 1 and p(r) = -— = p, (2.5) 

Po V 

with pQa^ = l/-\/2 + 3/(47r) being a rough estimate of the maximum density accessible to a 
disordered system [IS]. For a solid phase, the choice goes to 

o-l 3 



a 



1/3 



\PcpJ 



and p(r) = ^ 6'\r - Rj) , (2.6) 



where pcp is the number density at close packing {pcpO'^ = V^ for a FCC crystal) while 
{Rj} are the lattice sites. Upon making the further approximation of discarding the contri- 
butions to Fa from particles beyond the first coordination shell, one arrives at the following 



expression for the free energy per particle: 

keT [ln(pA3) - 1] - ksT In (1 - p/po) - 27re pa^ /+°^ dx x'^^{x) , for a fluid 

'^ ~ ^ A;bT [ln(pA3) - 1] - Sfc^Tln [l - [p/pcpf] - {zi/2)e<j) [{pcp/p)'/'] , for a solid 

(2.7) 

where A oc T~^/^ is the thermal wavelength and zi is the lattice coordination number. It 

is now straightforward to derive the entropy functions of the solid and fluid phases from 



Eq. (2.7). By eliminating T in favor of e, we obtain (up to an overall constant): 

kB ln(t> — vq) + {3/2)kB In (e + a/v) , for a fluid 

SkslnU/' - v'J^) + {3/2)kB\n\e+ {z,/2)e(l)\{v/vcpf' 



, for a solid 

(2.8) 

with a = 27iea^ J-^ dx x^(/)(x). 

In order to single out at a given T the most stable phase as a function of p, we should 



plot the two free energies in Eq. (2.7) vs. v = p~^ and then use the common-tangent 
construction. Alternatively, we may plot p as a function of P for flxed T and then look for 
(i) the crossing between the solid and fluid branches, and (ii) the "swallowtail" accompanying 
any isostructural phase transition if present. A typical outcome of this procedure can be 
seen in Fig. 1, showing the phase diagram for 0(x) = x^^ and zi = 12 (the solid is a FCC 
crystal). We see that a phase diagram of the standard simple-fluid type emerges in this case, 
which is enough for characterizing adiabatic freezing by the maximum-entropy method. 

Lastly, I describe the procedure by which the minimum of a convex multivariate function, 
here minus a total entropy, is computed. First a rough minimization of the objective function 
is attempted by the simulating-annealing algorithm [19] , which generates a random walk in 
state space which eventually brings to the sought minimum. For a convex function, this 
method is guaranteed to give the absolute minimum, that is the only minimum present, up 
to an error which decreases with the number of steps in the walk. Next, assuming that we 
got close to the extremum, a second optimization cycle is started with the gradient-descent 
method |20j, which eventually leads to the desired target state with high precision. 

III. RESULTS 

In the present Section, I report and carefully analyze the properties of the inhomoge- 
neous state attained by a supercooled liquid after its adiabatic relaxation to equilibrium, 



assuming the theory sketched in Section II. Two different experimental situations are dis- 
cussed, depending on whether a constraint is put on the system volume or pressure. In the 
constant-volume case, I consider the further possibility that a fixed amount of a foreign gas 
is present in the container. These cases are analyzed separately in the following. 

A. Constant volume 



Consider first a A^-particle liquid filling completely a closed rigid vessel of volume V. 
Initially at coexistence conditions, the liquid is subsequently driven metastable by slow 
cooling and then, after removal of the bath, violently perturbed in order to bring it to 
equilibrium. The final equilibrium state will maximize the total entropy ( 2.1[ ) (assuming 
no role for the vapor in the process, which is correct as long as the final pressure is larger 



than the triple-point value). Now, we specialize to a system described by Eq. (2.8), with 
(f){x) = x~^ and zi = 12. Initially, the hquid temperature and pressure are T^ and Pm, 
defining a point on the solid-liquid coexistence locus. In this state, the specific volume 
Vm = V/N can be obtained from 

ksTm 
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(3.1) 



When the liquid is brought to Tin < Tm at constant volume, its energy changes to 

_ 3 a 

2 Vm. 



(3.2) 



With these starting conditions, the state eventually reached by the system after disconnect- 
ing the bath and inducing solid nucleation is the one yielding the maximum of 
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(3.3) 



where ei = Ei/N,vi = Vi/N, and rii = Ni/N. The outcome of the maximization procedure 
are the thermodynamic variables characterizing the liquid fraction of the system in the final 
state. Any < n^ < 1 testifies of a partial crystallization of the liquid, hence of the stable 



coexistence of solid and liquid at equilibrium. Indeed, one easily derives from (2.1 ) that the 



necessary conditions for the maximum of Stot are the equahty of T, P, /i between the phases, 
as it may be checked a posteriori from the values of the temperature, the pressure, and the 
chemical potential of the solid and liquid fractions in the computed equilibrium state. 

I studied in detail the case Tm = 0.8 Tc. In Fig. 2 the final values (Tfin and Pgn) of 
the system temperature and pressure are reported as a function of Tin. Both quantities 
are smaller than the respective initial-state values, Tm and Pm', however, they too provide 
coordinates of points on the solid-liquid coexistence line. We also see that Tgn > Tin, i-e., 
the liquid heats up during the transformation, and the energy needed to the purpose clearly 
comes from the latent heat of freezing released during solidification. The solid fraction grows 
practically linearly with Tm — Tin (data not shown), up to about 20% for Tin = 0.1 Tm (note 
that there is no lower limit to undercooling in the present theory). As expected, the entropy 
increase in the transformation is larger the smaller Tin. 

B. Constant volume with a foreign gas in the vessel 

Let us now suppose that the liquid is prepared at Tm and Pm by exposure to a gaseous 
atmosphere (e.g. air), and that a small amount of gas gets trapped in the rigid vessel when 
sealing it. We then have a liquid in equilibrium with an immiscible gas in a container of fixed 
volume. For simplicity, I describe the gas as ideal and monoatomic, composed of Ng = XgN 



particles. The initial specific volume of the liquid, Vm, is still given by Eq. (3.1) but we now 
have Vm < V/N = t>tot- Specifically, 

-ftot = Vm + -^ ■ (3.4) 

As before, we imagine that the liquid and the gas are cooled very slowly until Tin is reached. 
At this point, the volume vi^ of the liquid is determined by minimizing the total Helmholtz 
free energy. This leads to the equation 



Vq Vf^ l^tot - ^'ir 



(3.5) 



which just represents the equality of pressures between the liquid and the gas. After removing 
the bath, we induce solid nucleation by a mechanical shock, and wait for the system to reach 



equilibrium. The state eventually attained is such as to maximize the total entropy 
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Upon maximizing Stot, one obtains the values of ei,vi,ni,es,Vs which provide a complete 



description of the equilibrium state. The five conditions for the maximum of (3.6) are 



equivalent to requiring the same temperature and pressure for the liquid, the solid, and 
the foreign gas in the final state (namely, Ti = Tg = Tg = Tgn and Pi = Pg = Pg = -Pfin); 
furthermore, also the chemical potentials of the liquid and solid fractions should be the same 
(i.e., fii = fig), indicating that the adiabatic decay of the metastable state eventually results 
in a stable coexistence between solid and liquid. 

In order to characterize adiabatic freezing, useful quantities to be monitored as a function 
of Tin are: the temperature and pressure of the solid-liquid mixture at equilibrium, Tgn and 
Pfin; the volume of the mixture, fmix = niVi + (1 —ni)vs, as compared to fin; and the entropy 
of the mixture, in comparison with the entropy of the supercooled liquid. For T^ = 0.8 Tc, 
I examined a number of Xg values in the range from 0.0001 to 1. For example, in Fig. 3 
the values of Tgn and Pfin are plotted for Xg = 0.001 and 0.1. Compared to the case where 
no gas is present, we see little differences for small to moderate undercoolings. However, 
below T ~ 0.45 Tm and as far as Xg <^ 1, we see a sharp change of slope in all curves, 
which is related to an abrupt crossover in the gas pressure at Tin (i-e., the right-hand side of 



Eq. (3.5)), from large to very small values. The crossover temperature T is roughly obtained 
by putting the left-hand side of Eq. (3.5) to zero for t>in = ftot (i-e., only below T, the gas 



volume at Tin is a significant portion of the total volume). It would be interesting to see 
whether a similar behavior is observed in a real liquid in the deeply supercooled regime. 
The pairs (Tffn, Pffn) for various Tin values and for Xg = 0.001 were reported on the phase 
diagram in Fig. 1, so as to confirm that the final equilibrium states are indeed coexistence 



states. In the top panel of Fig. 4, the solid fraction Ug = 1 — rii is plotted for Xg 



0.001 



and 0.1. It steadily increases with T^ — Tin, at an almost constant rate only provided Xg is 
not too small. In the panel below, the entropies of the mixture and the supercooled liquid 



(3.6) 



are reported. For the same two Xg values, Fig. 5 shows the final volume of the solid-liquid 
mixture, in comparison with the volume of the supercooled liquid. A clear crossing between 
the curves is found for a certain x^-dependent value Tx of Tjn. While above Tx the decay 
of the metastable state is accompanied with the system contraction, the opposite (i.e., an 
expansion) occurs below Tx . 

C. Constant pressure 

As a third example, I consider a supercooled liquid which relaxes to equilibrium under 
conditions that are simultaneously isobaric and adiabatic. This can be realized by conceiving 
a non-rigid and adiabatic boundary between the system and an environment characterized 
by the same pressure in all states (i.e., a volume reservoir). 

The energy function of the reservoir is clearly E^^Sr, Vr, Nr) = —PVr + f{Sr, Nr), where 
P is a constant and / is an unknown function. This is tantamount to say that the entropy 
function is of the form 

Sr{Er, Vr, Nr) = g{Er + PVr, Nr) , (3.8) 

for a convenient function g. When a system with entropy S{E, V, N) is in contact with a 
volume reservoir, the equilibrium state of the composite system is such as to maximize the 
total entropy 

5tot = SiE, V, N) + g{E,,, ~E + P{V,,, - V), Nr) . (3.9) 

It is easy to check from the latter equation that a necessary condition for equilibrium is 
that the system pressure be also P. If, moreover, the system boundary is adiabatic, the 
only way the system can exchange energy with the reservoir is pressure work, that is AE = 
—PAV. We now ask what is the total-entropy variation AS'tot resulting from the transition 
of the system of interest from an initial state, (i?in, V^m, A^), to a final state, (E'gn, Van, A^). 
Considering that i^gn + -PVfin = -Em + -PMn, we end up with 

A5tot = 5(Efin, l^fin, N) - S{E,^, Vin, N) = AS, (3.10) 

meaning that the final equilibrium state would also maximize the entropy increase of the 



system alone. The only residual variable in (3.10) is e.g. Van, while i^gn = -Em — -P(Van — H 



For a two-phase system in contact with a volume reservoir, the total entropy reads 
5tot = Si{Ei, Vi, Ni) + S,{E„ V, N-Ni)+ g{E,^, -Ei-E, + P{Vot -V-V), Nr) , (3.11) 
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prescribing the same pressure P for both phases at equihbrium. Assuming an initial state 
where only the liquid phase is present, and using the first law of thermodynamics to prove 
that Eifi^ + Esfi^ + PiVifi^ + Vsfin) = Eiiji + PVi^iri, the total-entropy increase is again reduced 
to the system-entropy increase, in turn given by 

AS = 5K^/,fin, ^,fin, Ni) + S,{Ei,,^ + PVi,,^ - E,,fi„ - Py,,fi, - PK,fi„, K,fin, N - Ni) 

- Si{Ei;,^,Vi;,^,N). (3.12) 



If we now write the necessary conditions for the maximum of (3.12), which eventually yield 
the unknowns -E/^fin, VJ^fin, A^;, V^^.tin, we find that they prescribe the same temperature and 
chemical potential for each phase (i.e., Ti = Tg and fii = fig), as well as the equality of both 
pressures with P {Pi = P^ = P). In explicit terms, the function to be maximized is 

>Jtot 
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where fin and Cin are the specific volume and energy of the supercooled liquid at Tjn, 

2 = Pm and Cin = -ksTm . (3.14) 

As before, the choice was made that T^ = 0.8 T^. At variance with the previous case 
where V was a fixed constant, an isobaric-adiabatic relaxation of the supercooled liquid 
to equilibrium now brings the system invariably to the original phase-diagram point, i.e., 
Tfin = Tm and Pfin = Pm, but in the modified form of a solid-liquid mixture, whose specific 
volume fmix is different from Vm and always larger than Vin (Fig. 6, upper panel). Once more, 
the solid fraction in the mixture is found to increase, to all practical purposes, linearly in 
Tm — Tin (Fig. 6, lower panel). I finally observe that the above results were successfully 
checked against the method for isenthalpic freezing described in [T]. 

IV. DISCUSSION 

Using a mean-field theory for illustrative purposes, I have shown that, under adiabatic 
conditions, a supercooled liquid transforms into a stable mixture of solid and liquid just 
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(3.13) 



for entropic reasons. However, until now thermal insulation was a mere hypothesis, and 
the question remains as to what conditions should be met in order that the decay of the 
metastable state can be treated as adiabatic also in the presence of the bath. Clearly, a two- 
stage freezing scenario could only be viable provided the alleged adiabatic step is guaranteed 
to conclude very quickly after the appearance of the first solid nucleus. Hence, there is no 
way a thorough analysis of adiabatic freezing can get around genuinely kinetic issues (rate of 
nucleation, growth velocity, etc.), which however lie outside the scope of a pure equilibrium 
theory (and of the present study as well). 

Possible hindrances to effective adiabaticity are of at least two kinds, one system-specific 
and another of a more general type. The inability to grow the solid phase rapidly would be 
typical of good glass formers, i.e., systems with sawtooth-like potential-surface topographies. 
If such a system is undercooled down to a temperature Tin which is only slightly above the 
glass-transition temperature, we expect that the adiabatic-freezing stage of solidification 
will be skipped altogether and conventional diathermal freezing (directed from the surface 
inward) will occur instead. A different and more basic form of kinetic bottleneck to adiabatic 
crystallization will be described below, after including in the description also the energy cost 
of the interface between the phases. 

For a two-phase equilibrium system which is both thermally and mechanically isolated 
from the environment, the total entropy can be written in the setting - originally devised 
by Gibbs himself ^T\ - where thermodynamic properties are attached also to the dividing 
surface a between the phases (see e.g. [221 123] )• The energy of a planar interface of area A 
can generally be written as Efj = TS^ + 'yA + fiNcr, where T and fi are those of the coexisting 
phases and 7 is the interface free energy (surface tension). By the Gibbs adsorption equation, 

S^dT + Ad7 + N^dfi = , (4.1) 

Ecr is reduced to just jA when the surface tension is independent of T and fi. Now switching 
to a spherical inclusion or droplet of the /3 phase in the metastable mother a phase, I 
make the further approximation that 7 is radius- independent, indeed a fair assumption only 
sufficiently close to coexistence (see e.g. [21]); moreover, I shall neglect surface-tension 
anisotropy, which is a small effect anyway for many crystals [22]. With these simplifications, 
the entropy of the a + P system becomes equal to: 

^tot(e;3, V(s, Np- E, V, N) = {N~ Np)s^{e^, i;„) + N^s^{ep, v^) , (4.2) 
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where ep and Vjs are the specific energy and volume of the nucleating phase, and Np is the 



number of particles in the droplet. In Eq. (4.2), the energy and volume of the mother phase 
are given hy Ea = E — Ej^ — E„ and Va = V — Vp respectively, or 

ii^, \^, A^ being the state variables of the composite system. For 7 = 0, the absolute maximum 



of (4.2)-(4.3) clearly coincides with the maximum of (3.3) for the given Tjn. 

The values of the internal variables ej3,Vj3,Np in a (possibly unstable) equilibrium state 
are obtained from equating the three partial derivatives of S'tot to zero. It is then a simple 
matter to show that these conditions are equivalent to: 

27 

/i/3(e/3,f/3) = fJ'a{ea,Va), (4.4) 

where r^ = (3A^/3f^/(47r))^/^ is the droplet radius. Hence, any cluster of the (3 phase which is 
in equilibrium with the a phase should have the same temperature and chemical potential as 
a, while the two pressures are different and related by the Laplace equation. In particular. 



Eqs. (4.4) would hold for the cluster of /3 phase in the inhomogeneous equilibrium state, 
associated with the absolute maximum of S'tot- 

Now take a to be the liquid (/) and (3 the solid (s), and assume these phases are de- 
scribed by Bans' theory. As far as the value of 7 is concerned, anything reasonable is 
good, for example the orientationally-averaged interfacial free energy of hard spheres 



7 = 0.561 ksTm/cr'^- I first checked that, upon maximizing (4.2)-(4.3) for a number of Tin 
values by the same numerical method as employed before, the conditions (4.4) are fulfilled. 
For small enough supersaturation, however, the maximum of the total entropy is invariably 
found at Ng = 0. In order to see what is going on, it is worth looking at the graph of 
the function AS = Stoti^s, Vg, N^; E, V, N) — Nsi{E/N, V/N), which represents the entropic 
advantage of the inhomogeneous system over the supercooled liquid. To simplify it further, 
AS* is projected onto the one- dimensional subspace where Cs and Vg are given the same 
values as in the point of absolute maximum of Stot- We are thus left with a function of A^^ 
only, which is reported in Fig. 7 for two small values of A^ (10^ and 10'') and a few under- 
cooling temperatures. A glance at Fig. 7 immediately reveals the existence of a sharp AS 
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maximum for a non-zero Ns value, corresponding to a two-phase equilibrium state. How- 
ever, a satellite maximum also exists at the origin, which is separated from the former one 
by an entropic "barrier" (the valley between the two peaks), and when the supersaturation 
becomes sufficiently small the absolute maximum of AS jumps to Ns = 0. Therefore, solid 
formation is thermally activated (i.e., it necessitates a favorable density fluctuation) and, 
for any fixed A^, there is a minimum undercooling threshold (which however is negligible 
for macroscopic A^) to overcome in order that solidification may occur. Below this thresh- 
old, the assumption of a rapid yet partial solidification of the liquid, which is at the heart 
of the present calculation, should be rejected - since no solid component, here modeled 
for simplicity as consisting of one single block, is found in the equilibrium state - and the 
onset of solidification is necessarily at the system surface in contact with the bath. Upon 
reducing the supersaturation further, the relative maximum for A^^ > disappears and no 
solid cluster, even only a metastable one, can form. A similar scenario is at work in the 
canonical-ensemble description of liquid nucleation from vapor [27] . 

Summarizing, the calculations in the present Section were aimed at checking whether 
the assumption of adiabaticity, which is at the basis of the results of Sections III, can 
survive the inclusion of the interface-energy contribution in the treatment. A necessary 
condition for that is a positive maximum of AS, which however only appears beyond a 
certain A^-dependent undercooling threshold, negligible in the large-size limit. This implies 
that small-sized liquids must be cooled sufficiently deep in order that freezing may start 
from the system interior; otherwise, homogeneous solid nucleation is obstructed (not simply 
activated!) and freezing will proceed diathermally from the outset, i.e., directly from the 
system boundaries. 

V. CONCLUSIONS 

Adiabatic freezing is the first lap of conventional freezing. It is observed whenever the 
energy released during solid nucleation in the very early stages of crystallization does not 
reach the thermostat but is almost completely spent in the heating up of the system, whose 
temperature raises quickly until solid-liquid coexistence is established at the equilibrium 
freezing temperature. Only later will crystallization proceed diathermally. 

I have studied the adiabatic freezing of a supercooled liquid using the van der Waals 
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theory of Refs. [T71 [18] , considering three possible experimental arrangements (constant vol- 
ume, constant volume with an entrapped foreign gas, and constant pressure). I have clarified 
that, when heat transfer to the external bath is kinetically hindered, the liquid undergoes 
partial crystallization just for thermodynamic reasons, i.e., as a result of total-entropy max- 
imization. Although the first outcome of recalescence is not usually the state of maximum 
entropy, due to the formation of a mushy zone which very slowly evolves to stable equi- 
librium, at least the trends exhibited by various system quantities with supercooling could 
roughly be predicted by simple thermodynamic arguments. 

In the attempt to unearth hidden hypotheses behind the modeling of the early stages of 
freezing as effectively adiabatic, I was finally led to consider the entropy of a liquid with a 
solid droplet inside. I have thus documented the existence, for a small-sized liquid system, 
of a minimum supersaturation to achieve in order that adiabatic freezing may occur under 
const ant- volume conditions. I defer to a future publication the real-life illustration of some 
of the features of adiabatic freezing that were highlighted in the present study. 
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FIG. 1: (Color online). Theoretical phase diagram for a system of particles interacting through 



the potential (2.2) with cl){x) = x~^ and zi = 12. Tc is the critical temperature and /3c = {ksTc)"^. 
The critical-point coordinates, pc and T^ follow from requiring that the first- and second-order 
density derivatives of the fluid pressure be simultaneously zero. One thus finds pc = Po/3 and 
ksTc = {8/27)apo, with a = (27r/3)e(7^. Top: phase diagram on the density-temperature plane, 
showing the extent of the coexistence regions; the triple temperature is between 0.6 and 0.65 of 
Tc. Bottom: phase diagram on the temperature-pressure plane, reporting as blue crosses also 
the (T, P) points characterizing the solid-liquid coexistence states borne out of the decay of the 
metastable- liquid states at various Tin values, for Xg = 0.001 (see Section III.B). 
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FIG. 2: Final equilibrium state after the adiabatic decay of the metastable liquid under constant- 
volume conditions, for Tm = 0.8 Tc. Top: temperature; bottom: pressure. 
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FIG. 3: Final equilibrium state after the adiabatic decay of the metastable liquid under constant- 
volume conditions, for T^ = 0.8 Tc and for two different amounts of foreign gas in the vessel 
(crosses, Xg = 0.001; squares, Xg = 0.1). Top: temperature; bottom: pressure. 
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FIG. 4: Top: Solid fraction in the equilibrium state resulting from the adiabatic decay of the 
metastable liquid under constant-volume conditions, for T^ = 0.8 Tc and for two different amounts 
of foreign gas in the vessel (crosses, Xg = 0.001; squares, Xg = 0.1). Bottom: Entropy of the 
solid- liquid mixture at Tgn (solid lines) vs. entropy of the supercooled liquid at T\a (dotted lines). 
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FIG. 5: Final equilibrium state after the adiabatic decay of the metastable liquid under constant- 
volume conditions, for Tm = 0.8 Tc and for two different amounts of foreign gas in the vessel (top 
panel, Xg = 0.001; bottom panel, Xg = 0.1). Volume of the solid-liquid mixture (solid lines) vs. 
volume of the supercooled liquid at Tin (dotted lines). 
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FIG. 6: Final equilibrium state after the adiabatic decay of the metastable liquid at constant 
pressure, for T^ = 0.8 Tc. Top: volume of the solid-liquid mixture at Tm (solid line) vs. volume of 
the liquid at Tjn (dotted line); bottom: solid fraction in the mixture. 
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FIG. 7: (Color online). Difference in specific entropy between the droplet-liquid mixture at Tgn 
and the original metastable liquid at Tin, as a function of the droplet "radius", A'^s . Two values 
of N are considered, 1000 (red curves, left) and 10000 (blue curves, right), for T^ = 0.8 Tc. For 
each N, various T\-a/Tc values were considered: from top to bottom, 0.57,0.60,0.63 for N = 10^; 
and 0.60, 0.65, 0.70 for N = 10^. 
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